High Temperature Superfluidity in Double Bilayer Graphene 
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Exciton bound states in solids between electrons and holes are predicted to form a superfluid at high tem- 
peratures. We show that by employing atomically thin crystals such as a pair of adjacent bilayer graphene 
sheets, equilibrium superfluidity of electron-hole pairs should be achievable for the first time. The transition 
temperatures are well above liquid helium temperatures. Because the sample parameters needed for the de- 
vice have already been attained in similar graphene devices, our work suggests a new route towards realizing 
high-temperature superfluidity in existing quality graphene samples. 
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It is proving a challenging task to observe superfluidity in 
a semiconductor electron-hole double quantum well system. 
Despite long-standing theoretical predictions Jl]|2] and signif- 
icant experimental efforts it is only very recently that 
Bose-Einstein condensation of excitons has been observed in 
this system. The transition temperature was low, ~ 1 K, 
and the condensate was non-equilibrium because of the fast 
recombination of the photo-excited excitons. Major obsta- 
cles blocking experimental realization of equilibrium super- 
fluidity in these systems include the following, (a) In most 
semiconductors, the electron and hole energy bands are badly 
mismatched. In GaAs not only do the effective masses differ 
by a factor of five, but the holes have highly non-parabolic 
energy dispersion and spin-| characteristics. [7j] (b) The large 
bandgap in GaAs means a bias of 1.5 eV needs to be applied 
across the dielectric barrier. This leads to a huge electric field 
across the barrier making it extremely hard to fabricate bar- 
riers thin enough for the electrons and holes to bind together 
while avoiding leakage between the layers. The formation of 
excitons is exponentially suppressed once the thickness Db 
of the barrier separating the two quantum wells exceeds the 
effective Bohr radius agQ]] (13 nm for GaAs). With the nar- 
rowest GaAs quantum wells widths Dw ~ 20 nm and barrier 
thicknesses Db ^ 10 nm to avoid leakage, the minimum ef- 
fective layer separation D e h — Db + Dw remains » aJ.lQl 

An alternative system, two graphene monolayers of elec- 
trons and holes separated by a dielectric barrier (2xMLG), has 
recently been proposed to observe this elusive superfluid. |@- 
rToh This system has some clear advantages over the GaAs 
system. Graphene is a gapless semiconductor with nearly 
identical conduction and valence bands so the mismatch be- 
tween the electron and hole Fermi surfaces is almost elimi- 
nated. In addition, the availability of very thin dielectric bar- 
riers separating the two monolayers makes the region with 
strong electron-hole pairing effects easily attainable. A barrier 
thickness as small as Db — 1 nm has already been demon- 
strated with a hexagonal boron nitride (hBN) dielectric. Ill 1 1 
The barrier can be made so thin for graphene both because 
there is no need for a large bias between the electron and hole 
layers, and also because hBN has a much larger bandgap (~ 5 
eV) than the barrier between double quantum wells in GaAs 



(~ 0.5 eV), allowing the barrier thickness to be readily re- 
duced without electrical leakage between the layers. 

However there is a new obstacle with graphene associated 
with the linear single-particle energy dispersion E±(k) = 
zLTiVFk. This makes it difficult to access the region of strong 
interactions. Also, bound excitons do not form because of the 
massless carriers. lll2ll The ratio r s = (Vcoui) / Ep is a use- 
ful measure of the importance of interactions relative to ki- 
netic energy. (Vcoui) = ie 2 1 K)^Trn is the average Coulomb 
interaction energy, where n is the charge carrier density and 
k is the dielectric constant. In graphene the Fermi velocity 
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is independent of n and the Fermi momen- 
tum for the spin-i carriers in 2D is kp = \j2irn/g v , where 
g v = 2 is the pseudospin factor. The Fermi energy is then 
Ep = Tivp y/2trn/g v . This makes r s = e 2 /(nhvp) < 1 for 
monolayer graphene, a constant of order unity, while recent 
theory suggests that an electron-hole superfluid can only oc- 
cur at measurable temperatures for r s > 2.3 in the 2xMLG 
system. II 13ll Very recent experiments have shown no evidence 
of superfluidity in this system despite achieving barrier thick- 
nesses as low as Db = 1 nm. ll ill This poses an exciting chal- 
lenge: can new experimentally realistic structures be designed 
using atomically thin crystals that allow the transition to a su- 
perfluid state? 

We concentrate here on bilayer graphene since it has been 
well characterized and exhibits extremely low levels of disor- 
der, but a number of other such crystals are possible. 11141 Our 
proposed system (2xBLG) consists of a pair of parallel bilayer 
graphene sheets (Fig.[TJ. The lower bilayer sheet is an electron 
bilayer comprising two parallel, A-B stacked, closely coupled 
electron layers of graphene with layer separation D e = 0.37 
nm. ll 1511 There is strong electron hopping between the two lay- 
ers in the sheet. The upper bilayer sheet is a hole bilayer 
consisting of two A-B stacked hole layers, but is otherwise 
analogous to the lower bilayer sheet with an identical layer 
separation Dh = D e . The two bilayer sheets are separated 
by a hBN insulating barrier of width Db to prevent tunneling 
between the sheets. There are separate electrical contacts to 
the two layers and a bias Vbb can be applied between them. 
The bias Vbb and biases Vtg and Vbg ° n top and bottom 
metal gates allow independent control over the carrier density 
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FIG. 1 : Spatially separated electron-hole system (2xBLG) with elec- 
trons in one graphene bilayer sheet separated by a hBN dielectric 
barrier from holes in a second graphene bilayer sheet. Top and bot- 
tom metal gates control the densities. 



in each layer and can adjust the symmetry of the electric field 
across the two sheets. By tuning the three biases, a wide range 
of electron and hole densities can be achieved. 

Bilayer graphene eliminates the problems caused by the lin- 
ear dispersion of monolayer graphene since over a wide range 
of electron or hole densities, 1 x 10 11 < n < 4 x 10 12 
cm~ 2 , symmetrically biased graphene bilayers behave as a 
zero gap semiconductor with a quadratic dispersion around 
the Fermi level|[ll: £±(k) ~ ±h 2 k 2 /2m*. The effective 
mass is m* ~ (0.03 to 0.05)m e , depending on the carrier 



density. [17] At very high densities n > 4xl0 12 cm 2 , E±(k) 



crosses over to linear behavior, but this density region lies in 
the weakly interacting regime and is not of interest. At very 
low densities n < n mm = 1 x 10 11 cm' 2 a trigonal warping 
of the bands transforms the two quadratic bands into three sets 
of Dirac-like linear bands. [ 1 8] This warping can be reduced if 
necessary by applying an asymmetric bias across each of the 
two bilayers which will also open up a gap separating the con- 
duction and valence bands. 

With quadratic energy dispersion, the Fermi energy in the 
bilayer sheet depends linearly on density, Ep = 7rh 2 n/g v m*, 
so the parameter r s now has a density dependence, r s = 
(e 2 / K)(g v m* /ft 2 )(l/y / 7rn). To experimentally reach large 
r s and the strongly interacting regime in the bilayer sheet, it 
is enough to decrease the carrier density n with suitable gate 
voltages. In graphene sheets the lowest density is restricted by 
onset of electron puddle formation, which occurs at densities 
below 10 11 cm -2 in high quality graphene bilayers. ill ill This 
corresponds to r s = 9 in the strongly interacting regime. 

We know that D e h -C a*, ensures strong electron-hole pair- 
ing. For a 2xBLG system embedded in a hBN dielectric with 
K = 3, lfl9ll setting m* = m* ~ mt = 0.04m e and using 
the reduced mass yields an effective Bohr radius a*, = 8 nm. 
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TABLE I: Parameters for 2xBLG, 2xMLG, and electron-hole GaAs 
double quantum well (GaAs DQW) systems. For a*, and Ry* we use 
the reduced mass. D e h is minimum effective layer separation, ro and 
r 3 are maximum particle spacing and maximum r 3 experimentally 
attained. (' Localized excitons do not form in 2xMLG. fl2lD 



This is large compared with the thickness of barriers akeady 
fabricated, Db — 1 nm ~ D„h. \\ ill 

Table I summarizes key physical parameters for three sys- 
tems in which superfluidity of spatially separated electron- 
hole pairs has been predicted. For GaAs DQW, although 
the E(k) is quadratic, the following experimental restrictions 
make it difficult to access the region with strong pairing ef- 



fects, (a) The barriers are wide, D e h ^ 2oq, 0, 12] (b) the 
largest value of r s attained is r s = 2, and (c) the effective Ry- 
dberg (Ry*) binding energy is small. For 2xMLG, the barriers 
can be very thin, D e h <C a*,, but the linear E± (k) keeps the 
system in the weakly interacting regime with r s < 1. What 
decisively favors the 2xBLG system is that the strong pair- 
ing regime is accessible in current samples because: (a) the 
extremely thin barriers Db > 1 nm and thin bilayer sheets 
D e = 0.37 nm, (b) the ability to tune r s to large values in or- 
der to access strong pairing, (c) the almost perfectly matched 
electron and hole bands resulting from their near equal effec- 
tive masses, leading to almost perfect particle hole symmetry 
and nesting between circular Fermi surfaces, and (d) the larger 
Ry* than for GaAs. 

We now calculate the superfluid energy gap and transition 
temperature to see whether a superfluid state can form in the 
2xBLG system with realistic sample parameters and at ex- 
perimental attainable temperatures. For simplicity we restrict 
Db > D e and approximate the system by a single layer of 
electrons I = e interacting with a single layer of holes t=h. 



The quadratic energy bands are £^ = 7ft 2 k 2 / (2m*) — //, 
where the band index 7=1 and — 1 for the upper and lower 
bands and // is the chemical potential. We make a standard 
transformation to positively charged particles for the holes. 
The effective Hamiltonian is 



fk 7 



k L k 



k-k' 



qk7 

k'7' 



7e f 7/1 f 7 e 7 h 



k+i 



~k+a-k' + ^ 



-k' + 



id) 



cZ and are creation and destruction operators for charge 
carriers in layer i and band 7. Spin indices are implicit. V^ k , 
is the static screened electron-hole interaction. 

A mean-field description is applicable on both the 
weak-coupling and strong-coupling sides of the BCS-BEC 
crossover for conventional pairing systems. [20] The mean- 
field equations at temperature T for // and the momentum- 
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FIG. 2: (a) Maximum of Ak at T = for different barrier thicknesses D e h- Above a critical density (arrows), A max drops discontinuously to 
sub-mK energies. Bands are not quadratic for n < n m i n . (b) Ak for D e h = 1 nm at densities (p), (q), and (r) marked in (a) 
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and /(-E) is the Fermi distribution function. Since we are in 
the strong coupling regime, we retain only the s-wave har- 
monic in the graphene geometrical form factor. [113 1 

We calculate the static screened V£ h (T) within the RPA, 



Kf(T): 



1 - 2v q Uo(q, T) + vffiiq, T)[l - e -^ D ^} 



(4) 



-2-ne 2 j{nq) is the unscreened Coulomb interaction. 



n (q,T) = (q, T) + n y u ' (q, T) is the sum of the normal 
and anomalous polarizabilities of a graphene bilayer sheet, de- 
fined in Eqs. 6 and 7 of the supplementary material. 

There has been considerable discussion on the effect of 
screening on superfluidity in electron-hole systems. Screen- 
ing weakens the interactions that drive the superfluid but there 
is an important subtlety: the appearance of even a small super- 
fluid gap in the excitation energy spectrum would completely 
suppress small momentum transfer screening. QlOl H3l I2JJ1 In 
the 2xMLG system, screening suppresses superfluidity for all 
physically accessible parameters. Ill 311 but in the 2xBLG sys- 
tem the interactions are stronger leading to more localized 
pairs and this can further reduce the effectiveness of screen- 
ing. Ultimately only full calculations can reveal whether the 
reduction is sufficiently great to allow onset of high-T super- 
fluidity, but the experimental observation in Ref. 6 of Bose- 
Einstein condensation of dilute electrons and holes in GaAs 
shows that screening does not invariably suppress superflu- 
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idity. Moreover, the fact that the measured transition tem- 
peratures in Ref. are consistent with mean-field calculations 
using unscreened Coulomb interactions, demonstrates that su- 
perfluidity can indeed suppress screening. 

We solved Eqs. [6] to [4] at T = 0. Figure |2a) shows 
A max , the maximum of Ak at T = 0, for different bar- 
rier thicknesses D e h- The electron and hole densities are set 
equal n e = = n so only the reduced mass enters the 
equations. I22II We take equal effective masses m* = 0.04m e , 
noting that a 25% difference between m* and m* h lfl7ll results 
in only a 10% change in the reduced mass and A max . We re- 
strict our density range n > n m in to ensure that Ef lies in 
the quadratic energy band. We neglect the small contributions 
from the negative branches of the bands. 

Figure |2ja) shows for each there is a critical density 
n c above which the gap A max is, at most, in the sub-mK en- 
ergy range. In realistic disordered systems it is unlikely there 
would be pairing in this case. At n = n c there is a sudden 
discontinuous jump in A max to much higher energies. Refer- 



131 reported a similar effect for the 2xMLG system but 



for r s > 2.35 only, and this cannot be achieved in the 2xMLG 
system. As D e h decreases, the pairing interactions become 
stronger and the superfluidity persists up to a higher n c . 

The nature of the superfluidity can be understood by look- 
ing at the k dependence of the gap Ak- In the weak-coupled 
BCS limit of Cooper pairs there would be a peak in Ak cen- 
tered at k = Jzf, while in the BEC region the peak in Ak 
would be centered at k = with a long tail falling off as 1/k, 
indicating pairs localized in real space. Figure |2jb) shows Ak 
at different densities for D e h = 1 nm. The densities (p), (q), 
and (r) are marked on Fig. |2ja). At density (p), Ak is con- 
stant out to k ~ 4:kp, indicating an average pair radius com- 
parable to the average particle spacing. Thus we are close to 
the BEC region. At densities (q) and (r) the peak is becoming 
less broad but is still centered on k = with a long 1/k tail. 
Further increasing the density should cause a peak at kp to 
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form but the gap collapses before that happens so the super- 
fluid does not reach the BCS limit before the effectiveness of 
the screening in Eq.|6]discontinuously jumps and superfluidity 
is suppressed. 

The very large gaps shown in Fig. EJa), in some cases 
A max > 500 K, are particularly interesting as, in contrast 
with high-Tc superconductors, these values of A max are main- 
tained across a wide density range. For reference, at n = 
5 x 10 11 cm~ 2 the Fermi temperature Tp = 175 K. 

To determine the superfluid transition temperatures T c , we 
first recall that for superfluids in 2D the mean-field criti- 
cal temperature T C MF , the temperature at which the mean- 
field superfluid gap goes to zero, can overestimate the true 
T c . The superfluid transition in 2D systems has a topological 
Kosterlitz-Thouless (KT) character|23i, and T c is determined 
by the KT temperature, 



T c = T KT = {ir/2)p s {T KT ) . 



(5) 



p s (T) is the superfluid stiffness. Mean-field theory gives a 
good estimate of p s (T) for quadratic bands in both the BCS 
and BEC limits.!!] p a (0) = E F /4ir at T = 0, and p s (T) falls 
off slowly up to T ~ A max if kpD e h is small, generally the 
case for our parameters. For n < n c we can take p s (T) ~ 
p s (0), Eq.[5]then giving Tkt = Ep/8. This is because even 
at a density as high as n — 6 x 10 11 cm -2 Ep/8 ~ 26 K, 
still much less than the A max shown in Fig. [2 a). For n > n c , 
A max is extremely small. Since p s (T) collapses to zero as T 
becomes larger than A max , T c = Tkt is similarly small for 
n > n c . 

Figure [3] shows the T-n phase diagram. For densities n 
greater than a critical density n c , the system is a Fermi liq- 
uid at all practicable non-zero temperatures. For n < n c and 
T below the transition temperature T c , the system is a super- 
fluid. T c = Tkt = Ep/8 grows linearly with density so 
the maximum transition temperature T" lax = Ep(n c )/8. For 
n < n c and T > T c there will still be strong signatures of 
the underlying superfluid state through the pseudogap. The 
pseudogap is a normal state precursor of the superconducting 
gap due to local dynamic pairing correlations and is produced 
by non-coherent fluctuations of the pairing field. Il24ll The de- 
viations from Fermi liquid behavior due to the pseudogap can 



persist up to room temperature, T ~ A max .|25|, 

Table HI] gives values of T™ ax and n c which increase with 
decreasing D e h because of the stronger electron-hole cou- 
pling. Using hBN barriers between the gates and between the 
graphene sheets (System (A)) gives k = 3, and T" lax is well 
above liquid helium temperatures. If the two graphene bilayer 
sheets are separated by a hBN dielectric barrier and suspended 
in air between upper and lower gates (System (B)), the effec- 
tive k = 1.5 and T™ ax is greatly increased. Devices similar to 
System (B) have recently been fabricated with a single bilayer 
graphene sheet suspended between two gates. 12711 To compare 
with existing literature Js, 13] on an idealized and unrealistic 
configuration of 2xMLG with k = 1 and D e h = 0, in 2xBLG 
the resulting T™ ax = 350 K, the maximum value of A max is 
7500 K, and n c corresponds to r s = 2.8. 
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FIG. 3: Density n - temperature T phase diagram (see text). 
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1 nm 
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TABLE II: Maximum transition temperatures T™^* and critical den- 
sities n c for barrier thicknesses D e h- System (A) is for the two 
graphene bilayer sheets embedded in a hBN dielectric. System (B) is 
for the two graphene bilayer sheets separated by a hBN barrier and 
suspended in air between the gates. 



We now discuss the experimental consequences of our 
phase diagram. Even at the lowest density n min = 10 11 
cm~ 2 , the transition temperatures in systems (A) and (B) re- 
main above liquid helium temperatures. There will be strong 
experimental signatures of the superfluid below T c . Although 
the pairs are neutral, the ability to make separate electrical 
contacts to the electron and hole layers allows for spectacu- 
lar electrical effects. [28] Coulomb drag and inter-layer tun- 
neling measurements will show significant enhancements as 
T is decreased below T c ,|0] while counterflow measurements 
can directly probe the superflow. [ 28j 29] Even up to room 
temperature there will be strong signatures of the pseudogap 
in phenomena such as compressibility, specific heat capacity 
and spin susceptibility. The wide density range over which 
superfluidity can be observed in a single device is in marked 
contrast with high-T c superconductors, where superconduc- 
tivity occurs only in a narrow 30% band of doping centered at 
optimal doping. 

The 2xBLG system is the first multiband system with just 
one condensate. It is an opposite case to multiband super- 
conductors such as magnesium-diboride where pairing is only 
within the bands and multiple coupled condensates appear, 
one condensate and one gap for each band. lT30ll In contrast, in 
the 2xBLG system, pairing is only possible between the dif- 
ferent hole and electron bands, leading to a single condensate 
and a single gap. 

The combination of extremely thin barriers, large r s , near- 
equal effective masses, and strong pairing attraction makes 
the 2xBLG system ideal for observing high-T c superfluidity. 
In existing quality samples it should be possible to study the 
phase transition from Fermi-liquid to superfluid as well as the 
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pseudogap physics. Our approach would also apply more gen- 
erally to other small gap semiconductors that can be made into 
atomically thin flakes and used with dielectrics such as hBN. 
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Methods 

We have taken the case of equal electron and hole den- 
sities (n e = n/j = n) and equal effective masses (m* = 
to* = nit). For zero temperature, the mean field Eqs. (3)-(4) 
in the main text for the momentum-dependent gap function 
= A^(T = 0) and chemical potentials pf and p h reduce 



to coupled equations for A^ and the average chemical poten- 
tial, fj,= (pi e + M/i)/2, 

n = 2 9v Y / (vr 1 ) 2 - (7) 

k 

£k = r y^i 2 k 2 /(2m*) — p, where m* is the reduced effec- 
tive mass of the electron-hole pairs. The screened zero-T 
electron-hole Coulomb attraction V£_ k , depends on the mo- 
mentum exchanged in the scattering process, thus making the 
gap momentum-dependent. The extra factor of i in Eq. (6) 

comes from the graphene geometrical form factor F^2, ■ In 
the phase space region of interest to us, coupling is strong and 
the s-wave harmonic contribution to F^, dominates, leading 
to the factor i. 

The numerical sums are transformed into continuous in- 
tegrals and then numerically evaluated using the Gauss- 
Legendre method. We numerically solved the T = non- 
linear system of integral equations, Eqs. (6) and (7), iteratively 
using a multivariable Newton method. The convergence ratio 
requirement for the variables was set at 10 -6 . Retaining only 
the s-wave harmonic of the form factor, results in a gap func- 
tion Ak that depends only on the modulus of the momentum 
k. 

Concerning the ultraviolet behavior of the mean field equa- 
tions, the integral over the momentum is well-defined, thanks 
to the exponential decay in momentum space of the bare 
Coulomb interaction VqC~ qD " h between electrons and holes 
in the bilayers separated by a dielectric barrier of thickness 
D e h- In our calculations we fixed the upper momentum cutoff 
for the integrations as k c = 5 x max(l/D e /j, Uf). Numerical 
stability was checked by increasing k c and increasing number 
of integration points. We confirmed that the results obtained 
by the Newton method for p > agreed with a solution of 
the mean field equations independently obtained by a direct 
recursive method. 



Suppression of screening: polarization bubble 

In the superfluid state at low temperature, the static limit 
of the polarization bubble in the diagrammatic (RPA-like) re- 
summation is suppressed at small momenta q. Reference 13 
in the main text discusses this effect in the limiting case of no 
dielectric barrier, D e h = 0, for energy bands with linear dis- 
persion. The polarization bubble is responsible for renormal- 
izing the bare Coulomb interaction. A small-q suppression of 
the screening allows strong unscreened electron-hole pairing 
peaked at small-q to occur, and this strong pairing can lead to 
a large gap in the excitation spectrum. 

We performed numerical calculations for the polarization 
bubbles for the case of quadratic energy bands. In the super- 
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FIG. 4: rig™' (q) and n o a ' (q) with two normal and two anomalous 
Green functions, respectively. 



fluid state the polarization bubble is given by the sum of the 
normal and anomalous bubbles (Eq. 4 in the main text and 
supplement Fig. 4). 

In RPA the polarization bubbles are constructed with the 
normal and anomalous Green's functions of BCS theory. As- 
suming the Green's functions are diagonal in the band indices, 
the normal and anomalous Green's functions are given by, 
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respectively, where u 7 , and v], are the Bogoliubov factors 



iulf 


= \ 






K 7 ) 2 


= \ 




a 

El 



(10) 



The summation of the fermionic Matsubara frequencies e„ of 



the particle-hole loop can then be done analytically, and taking 
the static limit the expressions for the bubbles are given by, 
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When the chemical potential /i is positive, the contributions 
of the lower band (7 or 7' = —1) to the total polarization 
bubbles are small because the corresponding fc-states are filled 
and far in energy from fi. Hence they contribute little to the 
particle-hole processes. We considered only the 7 = 7' = 
1 contribution and, as in the gap equation, only the s-wave 
harmonic of the geometrical form factor of graphene -FjJk ■ 

Physically, the effect of screening is suppressed because of 
the opening of the energy gap A at the Fermi surface. This ex- 
ponentially suppresses particle-hole processes with energies 
less than the gap energy A and it is precisely these low-energy 
processes that are needed to screen the long-range Coulomb 
interaction. Once a gap appears in the excitation spectrum, 
screening exactly vanishes in the q = limit. From a di- 
agrammatic point of view, in addition to a suppression of 
the bubble with diagonal Green functions at low momentum, 
there is an additional canceling contribution from the anoma- 
lous bubble with off-diagonal Green functions. At q = the 
cancelation is exact for any non-zero A. 



